Heatmap any DEG
library(pheatmap)
data = read.delim('~/DEGs_Infected_vs_Control_Combined.csv')
anycommondeg = data[data$Number_UP != 0 | data$Number_DOWN != 0,3:7]
anydeg = data[data$Number_UP >= 2 | data$Number_DOWN >= 2, 3:7]
rowd = c("correlation","euclidean", "maximum", "manhattan", "canberra", "binary", "minkowski")
cold = c("correlation","euclidean", "maximum", "manhattan", "canberra", "binary", "minkowski")
hclust = c("ward.D", "ward.D2", "single", "complete", "average", "mcquitty" , "median" , "centroid")
for (i in rowd){
for (k in hclust) {
pheatmap(as.matrix(anydeg[,2:5]),
col = topo.colors(100),
clustering_distance_rows = i,
clustering_distance_cols = i,
clustering_method = k,
main = paste0(i,k),
labels_col = c('Adults','Children','Mice - Feet','Mice - LN'))
}
}